Curvature condensation and bifurcation in an elastic shell 
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We study the formation and evolution of localized geometrical defects in an indented cylindri- 
cal elastic shell using a combination of experiment and numerical simulation. We find that as a 
symmetric localized indentation on a semi-cylindrical shell increases, there is a transition from a 
global mode of deformation to a localized one which leads to the condensation of curvature along 
a symmetric parabolic crease. This process introduces a soft mode in the system, converting a 
load-bearing structure into a hinged, kinematic mechanism. Further indentation leads to twinning 
wherein the parabolic crease bifurcates into two creases that move apart on either side of the line of 
symmetry. A qualitative theory captures the main features of the phenomena and leads to sharper 
questions about the nucleation of these defects. 

PACS numbers: 46.70.De,46.32.+x,46.70.Hg 
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The formation of defects in continuum physics has been 
the subject of long-standing investigation in many or- 
dered and partially ordered bulk condensed matter sys- 
tems such as crystals, liquid crystals and various quan- 
tum "super" -phases [lj. I n low-dimensional systems, the 
formation and evolution of defects involves an extra level 
of complexity due to the interplay between geometry and 
physics and remains an active area of research. Even in 
simple physical systems such as the mundane paper we 
write on and the textile sheets we wear daily, the vi- 
sual effects of these defects is arresting; in Fig^we see 
a heavy hanging drape with curved catenary-like wrin- 
kles that themselves have crease-like anisotropic defects 
(see arrows) on scales much smaller than the crease but 
much larger than the thickness of the drape. This hi- 
erarchy in structural scales is common in many elastic 
systems and naturally suggests the question of how such 
structures form. While the actual shape, response and 
stability of these structures has been the subject of much 
study recently 0], Q, 0, 0, following the early work 
of Pogorelov |6j , the process by which they arise from a 
featureless sheet is essentially unknown. 

To address this question in the simple context of an 
everyday example, a thin mylar sheet is bent into a half- 
cylindrical elastic shell (thickness t, radius R and length 
L, t <C R < L; R/t ~ 100) that is clamped along its 
lateral edges and indented at one end as shown in Fig. |3 
(a-d), a system on which preliminary investigations were 
conducted a while ago [7j. A point probe that is free to 
slip on the mylar is then used to indent the sheet at one 
end along the axis of symmetry. The shape of the sheet 
is then reconstructed using laser aided tomography [8j, 
wherein a laser sheet is used to interrogate the surface 
and determine its height. A CCD camera with a resolu- 
tion of 1024 x 768 pixels is used to image the reflected 
light and the laser sheet is rotated with a stepper motor 
which allows the entire surface to be scanned. 

For small values of the indentation 6, the shell deforms 




FIG. 1: Condensation of curvature in a heavy drape. The 
smooth catenary-like folds with a single direction of localized 
curvature coexist with regions of curvature condensation, in- 
dicated by arrows, where curvature is strongly focused. 



strongly in the immediate neighborhood of the edge, and 
the deformation decays monotonically on a scale com- 
parable to the cylinder radius (FigOi)- As the indenta- 
tion is increased further, we find that a localized zone 
of high curvature appears along the axis of symmetry, 
but at a distance d y away from the point of indentation 
(Fig|2h)- This "condensation" of curvature along a lo- 
calized parabolic defect in an elastic sheet is similar to 
that at the focus of a developable cone . For larger 
values of the indentation, if R/t is sufficiently large, the 
parabolic crease of localized curvature bifurcates into two 
such creases that move away symmetrically from the axis 
of the cylinder, in a manner analogous to defect twin- 
ning (Fig|5]:). Finally, for even larger values of the in- 
dentation, the defects stop moving along the axis of the 
cylinder, and a large region of the shell folds inwards 
(FigEJi). For long cylinders, the only intrinsic dimen- 
sionless parameter is the ratio R/t; there is effectively no 
dependence on material parameters since the long wave- 
length modes of deformation of a thin shell associated 
with bending and stretching are both characterized by 
the Young's modulus E and thickness t 0] . This makes 
the problem particularly attractive since geometry is at 
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FIG. 2: (Color online) Curvature condensation in an indented 
semi-cylinder. Figs (a) to (d) correspond to the experimen- 
tally determined deformation of the sheet with the colormap 
representing curvature h, yy along the axis of cylindrical sym- 
metry. Figs (e) to (h) show the evolution of the curvature 
h tVy with increasing indentation 8 (sheet viewed from the 
top) , for t/R — 0.001 determined using numerical simula- 
tions (see text); red corresponds to high curvatures, blue to 
low.. The Gauss curvature kg ~ \h yy shows the same qual- 
itative trends in both cases. 



the heart of all the observed phenomena. 

Understanding the experimental results requires the 
solution of the nonlinear equations for the deformations 
of shells 0. However, these are analytically insoluble 
except in some special cases involving axisymmetric ge- 
ometries and deformations. Therefore, we resort to a 
numerical investigation, and complement this with sim- 
ple scaling analyses to tease out the qualitative aspects 
of our results. A commercially available finite element 
code ABAQUS is used to minimizes the elastic energy of 
the shell (made of a material that is assumed to be lin- 
early elastic, with Young's modulus E = 100 MPa and 
Poisson's ratio v = 0.3) with an energy density 
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FIG. 3: Global and local modes of deformation, (a) Experi- 
mentally determined location of the condensate as a function 
of indentation depth, for a mylar sheet (t/R ~ .01). Filled 
symbols represent the condensate (before bifurcation) and 
open symbols represent the twinned state (after bifurcation), 
(b) We can understand this transition by probing the Gauss 
curvature along the axis of symmetry k,g(0,u) ~ ^h >yy (0,y) 
shown in determined numerically. We see the transition from 
a global mode with just a single curvature maximum (at the 
point of indentation) to a localized mode of deformation with 
a secondary maximum in the curvature that arises when the 
indentation S is larger than some threshold that is a function 
solely of the cylinder geometry. Here t/R = 0.001. 
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Here the first line characterizes the energy associated 
with in-plane deformations (with strains £1,62,7 ) and 
the second line accounts for the energy associated with 
out-of-plane deformations (k\,K2,t are the curvatures 
and twist relative to the undeformed tube). Four- node, 
quadrilateral shell elements with reduced integration and 
a large-deformation formulation were used in the calcu- 
lations. 

The results are shown in (Fig. EJe-h): as the indenta- 
tion is increased, the curvature condenses at a distance 
away from the edge, and this "condensate" bifurcates 
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into two upon further indentation. To understand these 
transitions qualitatively, we note that a complete cylin- 
drical shell that is pinched gently at one end responds by 
deforming globally in a nearly inextensional mode. Al- 
though a semi-cylinder clamped along its lateral edges 
also does the same when indented weakly, this global 
mode is screened by the lateral clamps. Strong inden- 
tation then leads to a transition from the global mode 
to a localized mode wherein the curvature localizes in 
a very small region at a distance d y from the edge. In 
FigEK? we show that this transition arises only when the 
indentation is larger than a threshold. To further inter- 
rogate this phenomenon, in FigOUs, we plot the Gauss 
curvature along the axis of symmetry of the cylinder as 
a function of the indentation. For small indentations, 
we find that the Gauss curvature along the axis of sym- 
metry decreases monotonically from the zone of inden- 
tation as shown in Fig^. However as the indentation 
is increased, the Gauss curvature develops an inflection 
point and eventually a secondary maximum, where the 
condensate forms. This condensation is crucially depen- 
dent on the two-dimensionality of the problem, and is 
not seen in one dimensional elastic systems |l 0f | . There 
are rough analogies to a phase transition when inter- 
preted in terms of the delocalized Gauss curvature and 
its localized condensate. Indeed the infinite-dimensional 
cylindrical structure with finite stiffness becomes a finite- 
dimensional mechanism, acquiring a new soft degree of 
freedom with a very low stiffness due to the formation of 
a hinge at the location of high curvature, which is akin 
to a condensed phase. 

As the cylindrical structure develops a hinge at loca- 
tion d y in response to the indentation (5, geometry implies 
that d y ~ 8. Of course as 8 is increased further, we ex- 
pect d y to eventually saturate at a location d VtSa t ~ R- 
In Fig. we see that numerical simulations confirm 
these naive expectations: for small indentations following 
the onset of localization, the location of the condensate 
increases linearly with 8 before it eventually saturates 
with dy ySa t ~ R- For very thin sheets, the curvature 
condensate bifurcates into two condensates when the in- 
dentation 8 > 5bif, each of the symmetric condensates 
moves away from the axis of symmetry, as can be seen 
in Fig. 0b. The critical indentation 8uf where the con- 
densate bifurcates into two is found to follow the scaling 

The location of the condensate is also a function of the 
cylinder thickness t; thicker, stiffer sheets lead to larger 
dy iSa t- In Fig^Jt, we see that d ViSa t is indeed larger for 
greater values of t since it is energetically favorable to 
have the curvature condense further from the indenta- 
tion. Zooming in on the region of curvature condensa- 
tion, we observe that the curvature is localized along a 
parabolic crescent-shaped structure (FigEt) similar to 
the core of a developable cone 0, 0, IE where stress 
is focused anisotropically. Our simulations show that for 
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FIG. 4: (Color Online) Curvature condensation and twinning, 
(a) The axial location d y (along the symmetry axis) of the 
condensate as a function of normalized indentation 8/R for 
different thicknesses (t/R as shown in the legend), with the 
dash-dotted lines showing saturation. The asterix represent 
experimental data from Q. The dashed line shows linear 
scaling in the initial stage, (b) The defect location d x as a 
function of the indentation follows the simple law d x ~ (8 — 
Sbif) 1 ^ 2 (dashed line) in the neighborhood of the bifurcation 
from one to two defects. 



a given indentation, the narrowest width of the parabola 
w ~ t occurs along the axis of symmetry, while its ra- 
dius of curvature R p ~ 2i 1 / 3 i? 2 / 3 (FigEb). In FigGJ we 
see that the scaled location of the curvature condensate 

d/R^ (t/R)- 1 / 3 . 

We now seek a qualitative understanding of these re- 
sults using the Doncll-Fdppl-von Karman equations that 
characterize the in and out -of-plane deformations of a 
weakly curved sheet In the neighborhood of the con- 
densate, the elastic bending energy Eg w Et 3 (^£j) 2 R P 2 , 
where 8 / dw is the dominant component of the curvature 
of the shell along the cylinder axis, and R p is the approx- 
imate area of inversion of the cylinder due to the forma- 
tion of the crease of radius R p . Stretching in the vicinity 
of the crease is primarily due to the non- vanishing Gauss 
curvature of the crease kg ~ -^^j via Gauss's Theorema 
Egregium. Then, we may write the stretching energy 
Es ~ Et-f 2 R p w , where the stretching strain 7 ~ cr/Et 
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FIG. 5: (Color Online) Length scales characterizing the con- 
densate as a function of the cylinder geometry t/R. (a) A 
magnified view of the curvature condensate in Fig. |5f show- 
ing the lengthscales R v and d. (b) The normalized radius of 
curvature R p /R of the parabolic crease as a function of the 
normalized thickness t/R, at 5 = 0.08, with the solid line 
~ (t/R) 1 / 3 . (c) The normalized location of the condensate 
d/R as a function of t/R at 8/R = 0.05, with the solid line 
~ (t/R)~ 1/[i . In each case, the open circles correspond to 
numerical simulations. 

in terms of the in-plane stress a ~ E£\J~ 2 kq and R p w 
is the area over which this stretching arises. This last 
relation follows from the compatibility relation for the 
stresses in the Foppl-von Karman equation which reads 
V 2 er ~ —kg- The narrowest width of the anisotropic 
crease is limited by the thickness so that w ~ t, consistent 
with our numerical simulations. Using this fact in min- 
imizing the total energy Eb + E$ yields R p ~ t 1 / 3 i? 2 / 3 
which agrees with our simulations over the range explored 
(FigOa) , consistent with earlier results on the size of the 
core of the defect in a sheet bent into a cone 0, 0, |n| . 
Assuming that indentation leads to an effective deforma- 
tion wherein the cylinder is reflected about a parabolic 
boundary in the neighborhood of the localized defect, 
we find that d ra 8R/R p ~ S(R/t) 1 ^ 3 , seen in the simu- 
lations (FigEJ;). For very thin shells, the curvature con- 
densate bifurcates into two which move to a location ±d x 
on either side of the axis of symmetry as the indentation 
increases beyond a threshold. Indeed, symmetry con- 
siderations thus demand that this supercritical pitchfork 



bifurcation must follow the scaling d x ~ (5 — S^f) 1 / 2 con- 
sistent with what we observe in our numerical simulations 
(Fig. Eb). 

Our studies here are but a first step in understanding 
the question of defect formation and evolution in an elas- 
tic system. By indenting the edge of a semi-cylindrical 
shell we followed the exploration of the geometric local- 
ization of deformation experimentally and numerically. 
Although we can understand some aspects of the curva- 
ture condensate, much still remains to be done in try- 
ing to quantify this phenomenon which is fairly generic 
in curved shells. For example, our understanding of 
the condensation of curvature is qualitative - we do not 
know the conditions under which this would happen more 
generally. In addition, our description of twinning re- 
mains qualitative being based on symmetry arguments, 
not mechanistic ones. Finally our studies suggest pos- 
sible applications to nanotechnology (e.g. deformation 
of nanotubes) , molecular biology (e.g. virus shells), and 
physical chemistry (e.g. morphology of colloidal parti- 
cles); indeed wherever these thin- walled structures arise. 
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